Unexpected non-Wigner behavior in level-spacing distributions of 
next-nearest-neighbor coupled XXZ spin chains 
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The level-spacing distributions of XXZ spin chains with next-nearest-neighbor couplings are stud- 
ied under periodic boundary conditions. We confirm that integrable XXZ spin chains mostly have 
the Poisson distribution as expected. On the contrary, the level-spacing distributions of next-nearest- 
neighbor coupled XXZ chains are given by non-Wigner distributions. It is against the expectations, 
since the models are nonintegrable. 
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I. INTRODUCTION 

Random matrix theories have been successfully applied 
to the energy spectra of not only chaotic systems but also 



quantum spin systems 



1.2.3.4.5.6.7.8.9 



If a given Hamilto- 



nian is integrable by the Bethe ansatz, the level-spacing 
distribution should be described by the Poisson distribu- 
tion: 



^'poi(s) = exp(-s). 



If it is not integrable, the level-spacing distribution 
should be given by the Wigner distribution: 



TTS / ITS 
Avig(s) = yCXp I — 



These behaviors of the level-spacing distribution have 
been observed in one-dimensional (lD))2ii& 2T)^i^ and 
SD*" spin systems. They are also confirmed for strongly 
correlated systemSfiSi and applied to the recent study of 
quantum dotsiii 

The numerical observationsiSii^iLSiifliii are impor- 
tant. In fact, for the quantum spin systems, there has 
been no theoretical or analytical derivation of the sug- 
gested behaviors of the level-spacing distribution. Ran- 
dom matrix theories do not necessarily extend to them. 
However, it seems that no counterexample has been 
shown explicitly for quantum spin systems. The sug- 
gested behaviors have been numerically confirmed for 
quite a large number of quantum spin systems and sta- 
tistical lattice systems iSiSi^iLSiS Moreover, the behaviors 
are often considered as empirical rules, by which we can 
practically determine whether a given lattice system is 
integrable or no^. Thus, it could be nontrivial if one 
would find such a quantum spin system that does not 
show the standard behaviors of the level-spacing distri- 
bution. 

Recently, it has been found that an extraordinary sym- 
metry appears for special cases of the integrable XXZ 



spin chain: the XXZ Hamiltonian commutes with the 
sh loop algebra at some particular values of the XXZ 
coupling associated with roots of unity>iS The loop al- 
gebra is an infinite-dimensional Lie algebra, which holds 
only at the particular values of the XXZ coupling and 
not for generic values. It is suggested numerically-'^'^ that 
the standard Bethe ansatz does not hold for the spe- 
cial cases. Furthermore, the dimensions of degenerate 
eigenspaces of the loop algebra are given by some expo- 
(X) nential functions of the system sizep-** and they can be 
extremely large. It should therefore be nontrivial how 
the large degeneracies are resolved by nonintegrable per- 
turbative terms in the spin Hamiltonian. Thus, the SI2 
loop algebra symmetry may motivate us to reconsider not 
only the Bethe ansatz solvability of the integrable XXZ 
chain but also the standard statistical behaviors of en- 
ergy levels for various XXZ chains close to the integrable 
one. We note that some level crossings of the sh loop 
algebra are shown in Fig. 1 of Ref. 



(2) 



In this paper, we discuss level-spacing distributions for 
the spin i XXZ chains on finite sites under periodic 
boundary conditions. We mainly discuss nonintegrable 
cases, although the present study has been originally mo- 
tivated by the s/2 loop algebra symmetry of the integrable 
XXZ chain. We first confirm that integrable spin chains 
show the Poisson distribution for generic cases. Here 
we exclude the special cases related to roots of unity. 
Then, we consider the next-nearest-neighbor (NNN) cou- 
pled Heisenberg spin chain (or XXX spin chain), which 
is nonintegrable. We now note that the XXX chain with 
NNN couplings has the Wigner distribution, as shown in 
Ref. 13- For the XXZ spin chains with NNN couplings, 
however, we find that the level-spacing distributions are 
not given by the standard Wigner distribution. The ob- 
servation should be nontrivial since the systems are non- 
integrable. Finally, we discuss possible reasons why the 
non-Wigner distribution is obtained for the NNN coupled 
XXZ spin chains. 
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II. NUMERICAL PROCEDURE 



The Haniiltonian matrices can be separated into some 
sectors; in each sector, the eigenstates have the same 
quantum numbers. This is because the eigenvalues with 
different symmetries are uncorrelated. The XXZ chains 
have three trivial symmetries, namely, reflections, trans- 
lations, and spin rotations around the z axis. Therefore, 
to desymmetrize the Hamiltonians, we consider three 
quantum numbers: parities, momenta, and total S^. 
We calculate the eigenvalues of the largest sectors. The 
largest sector of each Hamiltonian has 440 eigenvalues for 
the lattice size L — 16. To calculate the eigenvalues, we 
use standard numerical methods, which are contained in 
the LAPACK library. 

To find universal statistical properties of the Hamilto- 
nians, one has to deal with unfolded eigenvalues instead 
of raw eigenvalues. The unfolded eigenvalues are renor- 
malized values, whose local density of states is equal to 
unity everywhere in the spectrum. In this paper, the 
unfolded eigenvalues obtained from the raw eigen- 

values Ei in the following method. Let us define the 
integrated density of states as 



N 



n{E) ^Y.e{E- E, 



(3) 



1=1 



Here 0{E) is the step function and N is the number of 
the eigenvalues. We choose some points of coordinates: 
{Ei,n{Ei)) for i = 1, 21, 41, • • • , A^. The average of inte- 
grated density of states {n{E)) is approximated by the 
spline interpolation through the chosen points. The un- 
folded eigenvalues are defined as 



Xi = {ii{Ei)). 



(4) 



The level-spacing distributions are given by the proba- 
bility function P{s), where s — — Xi. 



III. NEXT-NEAREST-NEIGHBOR COUPLED 
XXZ CHAIN 

Let us introduce the spin i XXZ chain with NNN cou- 
plings on L sites by 



H = Ji ^ [SfSf^, + SfSl, + Ai5f 5f+,) (5) 
1=1 

L 



1=1 



where = (l/2)cr" and (cr"^, cr^, cr^) are the Pauli ma- 
trices; periodic boundary conditions are imposed. For 
simplicity, we put Ji = 1 hereafter in the paper. The 
Hamiltonian is nonintegrable for J2 7^ 0, while it is inte- 
grable for J2 = 0. 
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FIG. 1: Level-spacing distribution of the integrable XXZ 
chain (J2 — 0) for L — 16, Ji — 1, and Ai = 0.1 . Broken 
line shows the Poisson distribution. There is no degeneracy 
at s = 0. 



Let us confirm the Poissonian behavior for the generic 
case of the integrable XXZ spin chain. When J2 = 0, the 
level-spacing distribution P{s) mostly shows the Poisson 
distribution as shown in Fig. ^ We confirmed numeri- 
cally the standard result for some generic values of the 
XXZ coupling Ai (0 < Ai < 1). We note that we ex- 
clude the special values of the XXZ coupling that are 
given by Ai = cos(TO7r/A^) for some integers m and TV, 
where Ai is related to a root of unity, q, through the 
relation Ai = {q + l/q)/2}^ 

In nonintegrable cases (J2 7^ 0), however, P{s) shows a 
non-Wigner distribution against expectations. The non- 
Wigner behavior of the NNN coupled XXZ spin chains 
is generic. We have calculated P{s) for various combi- 
nations of J2 and A2 when a — 1 and for various val- 
ues of A2 when a = 0. Here, Ai = 1 for simplicity. 
When a — 1 and A2 = 0.98 , P{s) looks similar to 
the Poisson distribution rather than the Wigner distri- 
bution for any value of J2, as shown in Figs. Ela) and 
|2Ic). The Hamiltonian is close to but not exactly the 
same as the Heisenberg chain, when A2 = 0.98. We find 
such a Poisson-like distribution when A2 is very close to 
1.0 such as more than about 0.95. For the NNN coupled 
Heisenberg chain (A2 = 1), however, it was shown that 
P{s) has the Wigner distributioniSi^ 

When A2 decreases, P{s) becomes a non-Poisson dis- 
tribution as shown in Figs.HJb) and|3d). However, it is 
not a Wigner-like distribution, either. The level-spacing 
distribution P{s) strongly depends on A2. On the other 
hand, it seems that P{s) is almost independent of J2. 

When only one parameter is nonzero, namely, for the 
case of either a = or A2 = 0, the nonzero parameter 
does not change P{s) very much. In that case, P{s) is 
given by neither a Poisson-like nor Wigner-like distribu- 
tion. 

For all the investigated nonintegrable cases except for 
the case where A2 is close to 1.0 and J2 7^ 0, P{s) can be 
approximated by the arithmetic average of the Poisson 
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FIG. 2: Level-spacing distributions of NNN coupled chain for L = 16, Ji = Ai = q = 1, (a) J2 = 1, A2 = 0.98; (b) J2 = 1, 
A2 — 0.5; (c) J2 = 2, A2 = 0.98; (d) J2 — 2, A2 — 0.5. Broken lines, the Poisson distribution; dotted lines, the Wigner 
distribution. There is no degeneracy: P{s) = at s = for all the four cases. 



and Wigner distributions: 



IV. DISCUSSION 



Pav(s)= [Ppoi(s)+Avig(s)]/2, 



(6) 



as shown in Fig. |3| 



Let us remark on the special case of the integrable 
XXZ chain related to a root of unity. At the special 
values of Ai, we find that P(s) shows a novel peak at 
s = as shown in Fig. 01 For L = 16, large degen- 
eracies remain for Ai = cos(7r/i), i = 2,3,4, even after 
the desymmetrization procedure mentioned in Sec. ^1 is 
performed. We recall that for generic values of Ai, no de- 
generacy remains after the desymmetrization procedure 
is completed and the Poisson distribution is obtained as 
shown in Fig. The appearance of the peak in P{s) 
at s = should be consistent with the SI2 loop algebra 
symmetry, which holds for the transfer matrices of the 
XXZ and XYZ spin chains only at the special values of 
AiiiSiiii^ The result could be related to the observation 
in Ref. that the eigenvalue spacing distribution P{s) 
has a peak at s = for the transfer matrix of the eight- 
vertex model under the free-fermion conditions. Details 
should be discussed in forthcoming papers. 



Let us discuss some possible reasons why the noin- 
tegrable models give the non- Wigner distributions. We 
may consider two reasons: extra symmetries or finite-size 
effects. 

The non- Wigner distributions, such as shown in Figs. [21 
and 121 suggest that the Hamiltonian should have some 
extra symmetries. It must be a possible standard in- 
terpretation according to Refs. and 0. However, for 
the NNN coupled XXZ chains, it is not clear whether 
there indeed exist some extra symmetries other than re- 
flections, translations, and spin rotations around the z 
axis. 

We may also consider some finite-size effects, since the 
distribution P(s) is always obtained for some finite sys- 
tems. All the calculations in this paper are performed 
on 16-site chains. Recall that the largest sector of the 
NNN coupled XXZ Hamiltonian has 440 eigenvalues. 
The number is not small. If the level-spacing distribu- 
tion calculated on hundreds-site or thousands-site chains 
could give a Wigner-like distribution, then the finite-size 
effects should be nontrivial. 
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FIG. 4: Level-spacing distribution of integrable XXZ chain 
(J2 = 0) for L = 16, Ji = 1, Ai = cos(7r/3) = 0.5. Broken 
line shows the Poisson distribution. The peak at s = is 
given by degenerate eigenvalues. 



FIG. 3: Level-spacing distributions of NNN coupled chain 
for L = 16, Ji = Ai = J2 = 1, (a) a = 0, A2 = 1; (b) a = 1, 
A2 = 0. Broken lines, the Poisson distribution; dotted lines, 
the Wigner distribution; long dashed lines, the arithmetic av- 
erage of the Poisson and Wigner distributions. There is no 
degeneracy for both (a) and (b) [P{s) = at s = 0]. 
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